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Abstract 

I present a modification of the shadow-lattice technique, which allows one to 
derive low temperature series for discrete spin models to high orders. Results are 
given for the 3-d Ising model up to 64 excited bonds, for the 4-d Ising model up to 
96 excited bonds and the 3-d Potts model up to 56 excited bonds. 



Series expansions are a major tool to obtain critical properties of statistical systems and to 
test Monte Carlo programs. Usually series are derived using graphical techniques, but recently 
a method using direct enumeration of the relevant configurations on the lattice has been pro- 
posed Whereas the diagrammatic techniques have the disadvantage that they are hard 
to port to parallel or vector machines, the direct enumeration relies heavily on fast machines, 
taking full advantage of parallel architecture, but is computationally too demanding to be 
implemented on smaller machines. One advantage of direct enumeration is that it is done au- 
tomatically, i.e. once the program is set up it can in principle calculate the series to any given 
order. 

In this paper 1 shall sketch a method which is a variant of the shadow method (see for 
example 0) and which seems to be powerful, especially for higher dimensional systems. A 
detailed description will be published elsewhere 0. Although it is a graphical technique, the 
generation of the series is fully automatic, minimizing possible errors due to missing diagrams. 
Using this method, existing series for the Ising and Potts model @] [§] are extended by more 
than ten terms using 1 day of CPU-time on a personal computer. To obtain the magnetisation 
and free energy for the d = 3 Ising model up to 48 excited bonds 1 need 43 seconds, which has 
to be compared to more than a day on a Cray-YMP, which was needed in to obtain the 
free energy to the same order. Because the method uses an embedding which corresponds to 
real configurations, the calculation of correlation functions could be added rather easily and, 
based on the experience of 0, 1 expect an additional factor of 15 for the CPU-time needed. 
The method could be reasonably well adapted to vector computers, allowing the generation of 
series for various discrete models with a small amount of CPU-time. 

In the next section 1 shall describe the proposed method using the Ising model as an example. 
In section ^ the generalization to other discrete models is addressed. Finally some results are 
presented. 
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Consider the Ising model on a (i-dimensional hypercubic lattice of volume A'" in the presence of 
an external magnetic field h. The spin variables on each site a{x) can take on the values ±1 
and the action of the model is given by 

Sh = -2k J2 (^{^)(^{y) - hJ2^i^)^ (1) 

<xy> X 

where the first sum runs over nearest neighbours. 
The free energy per unit volume / is given by 

= - hm i-lnZ,. (2) 

iV— >(X) iV 

At couplings K ^ Kc most of the spins are aligned even with no magnetic field. Thus one 
can start with a lattice of N aligned spins and treat excitations arising from groups of fiipped 
spins. For zero magnetic field all states are doubly degenerate due to the symmetry a — > —a. 
One can overcome this degeneracy by introducing a small magnetic field. With u — e~^'^ and 
y — e~^'* the partition function for N spins can be rewritten as 

N qN ^ 

IuZn = -—Iny - —\nu + 22Pr,N{y)u'\ (3) 

^ ^ r->0 

where the Pr,N{y)u^ are contributions from fiipped spins. 

The task of the low temperature expansion is to obtain the Pr,N{y), giving an expansion of 
thermodynamic quantities in u. Because the free energy is an extensive quantity it suffices to 
calculate 

f = Epriyy^ (4) 

r>0 

where Priv) is the part of Pr^^iy) linear in N. 

One can divide the lattice into an even and an odd sub-lattice, depending whether the sum 
of coordinates for a given site is even or odd. The shadow method uses the fact that spins 
on odd sites do not interact with each other: If one flips a spin on an odd site, the change 
in the power of u is given by the number of its flipped neighbours, which are entirely on even 
sites. It is therefore sufficient to embed sub-configurations consisting only of spins on the even 
sub-lattice and note the number uj{o) of flipped neighbours for each odd site o. From now on 
I shall refer to odd sites as "shadows" and use the term "spin" exclusively for even sites. We 
can break up the free energy into contributions fs,t coming from s flipped spins and t flipped 
shadows: 

/ = E/m- (5) 

s,t 

Because the even and the odd sub-lattice are symmetric, one must get the same result when 
interchanging flipped spins and flipped shadows: 

fs,t — ft,s- (6) 

This property is referred to as code-balance. The advantage of the shadow method is that one 
has to embed approximately only half the spins with httle overhead for the bookkeeping of the 
uj{o). In practice one saves less CPU-time than expected, because usually all configurations up 
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to a given number s of spins are generated. The resulting series are correct up to the lowest 
order one obtains with s + 1 spins. Most of the s-spin contributions have a degree higher than 
that and are thus generated in vain. To give an example: only 1 out of 113636 (or even more) 
8-spin diagrams is needed to calculate the d = 4 Ising model series up to order 32. 

Here I propose to use code-balance in conjunction with predetermining the lowest degree 
of a diagram in order to generate only those diagrams which are actually needed for a given 
maximum order Imax in u. Because code balance is employed one needs to find all diagrams M 
which will give non-zero contributions fs,t{M) with s < t and an order in u less than or equal 
to Imax- Contributions fs,t with s > t are given by code-balance. 

Let me briefly review the steps necessary to construct a series: A configuration on a lattice 
will consist of one or more clusters of flipped spins where the different clusters are not joined 
by nearest-neighbour bonds, but the spins within a cluster are connected among themselves. 
The configurations are described by connectivity matrices. Because the shadow method is 
used, only the spins on even sites are entered into the connectivity matrix. An entry in the 
connectivity matrix describes whether two spins have 0, 1 or 2 shadows in common. Let me call 
a matrix describing s spins an s-matrix. It is impractical to embed a configuration consisting 
of more than one cluster because there will be of the order A^*^ different relative positions for c 
clusters. Instead one constructs only connectivity matrices which correspond to a single cluster 
and determines in a second step the number of ways two or more of such matrices can be put 
on the lattice. 

The construction of the connectivity matrices can be done recursively: One starts with the 
(trivial) matrix for one spin and adds another spin, which must be connected to the first spin 
by a non-zero entry in the connectivity matrix. Other spins are then added in such a way that 
each new spin is connected to at least one spin of the old matrix. 

Having constructed a connectivity matrix, one can embed it on the lattice in order to find 
its embedding number and its code (that is the set {uj{o)} for all shadows). 

The embedding number of two clusters on the lattice is the product of the embedding 
numbers of the two clusters (modulo symmetry factors) minus the number of embeddings 
where the two clusters are positioned in such a way that they actually form one cluster. Only 
the second term is of interest because the first term is ~ A^^ and does therefore not contribute 
to the free energy. Instead of trying to figure out when two clusters overlap, one can work the 
way backwards and determine how often a single cluster can be split into two or more clusters. 

As said above, one needs to consider only configurations with s < t. Consequently I define 
the degree of an s-matrix M as the minimum order in u of its contributions fs^t{M) with t > s. 
The calculation of the minimum degree can be performed easily if the code of the matrix is 
known. 

Because all the matrices considered are finite one can show that in each matrix there is at 
least one spin which can be removed without increasing the degree of the matrix and without 
disconnecting the matrix into more than one cluster. Turning this around, all required {s + 1)- 
matrices can be constructed by adding a spinQ to s-matrices, and the (s + l)-matrices will 
have at least the degree of the s-matrix they were constructed from. This allows us to discard 
immediately any s-matrices which have a degree higher than Imax- This in turn will prevent 
many (s + l)-diagrams of degree > Imax from being generated and will save a large amount of 
computational effort. For some s = Smax no matrices with degree less than or equal to Imax 
are produced when adding another spin. At this point all connectivity matrices have been 
constructed. 

^The conditions this spin has to fulfiU will be given elsewhere. 
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As was mentioned earlier I take a single cluster apart in order to obtain disconnected 
diagrams. Also for this procedure one can show that the degree of the resulting diagram is at 
least the same as the degree of the diagram one uses as starting point. 

Using these two observations one is able to construct the required series without generating 
unnecessary matrices or partitionings. 

Despite using code balance to construct the series, it can still serve to check the calculation: 
A matrix which contributes to the order of fs,t with s > t is guaranteed to be among the 
constructed matrices if 

k<lmax-{s-t){d-l). (7) 

This allows one to check ft^g against fs^t for orders k bounded by eq. (0). 

I performed the code balance test on series where some diagrams or some partitionings were 
deliberately missing. In all cases the series did not fulfill the code balance, giving confidence in 
the reliability of this test. But, as obvious from eq. (^, the highest {d — 1) orders of a series 
cannot be checked. 

3 

As an example of another discrete model let me consider the 3-state Potts model on a d- 
dimensional lattice with sites. The action reads: 

sh = -^Y. m^),v{y))-h{Y.m^)A)-\N] (8) 

<xy> \ X / 

with p{x) e {1,2,3}. 

Again one starts with a ground state where all spins are aligned with the external magnetic 
field. Using u = exp(— k), y = exp(— /i), the free energy for N spins can be rewritten as 

f = jN\ny + ^\nu - E^n^^'^*, (9) 

where the Pi^r,N are configurations with / left-pointing {p{x) = 2) spins and r right-pointing 
{p{x) = 3) spins. 

The expansion is done similarly to the Ising model. One only has to insert another step 
which selects a subset of the flipped spins to be left-pointing. The shadow method can be used 
as well, but instead of noting for each shadow the number of flipped neighbours, one has to 
note Ui, the number of left-pointing and Ur, the number of right-pointing neighbours. 

The minimum degree for a given s-matrix can be achieved by aligning all flipped spins, for 
example letting them point to the left. In this case the degree is equal to the one for the Ising 
mode^. Therefore the set of matrices is identical to the one for the Ising model. As in the 
Ising model, the series can be checked by code balance. 

4 

The number of matrices contributing for a given degree of the series are given in table |l]. The 
time needed on a personal computer doing about 12 MIPS is given in the last row of the table. 
Memory requirements were below 3Mbyte. For the Potts model the time needed was 300 sec for 

^Up to a factor which comes from the different definitions of u. 
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^max = 48 and 6775 sec for l^ax = 56. This is about a factor 6 more than for the corresponding 
Ising modeL 

Although the time needed for the d=3 Ising model is of the order of a day, a further extension 
of the series seems rather hard with this method. Even if one assumes a 1000 MIPS machine, 
I would expect 1 day of CPU time to go up to u^^ and a year for 

The expansion gives the free energy and all susceptibilities at the same time, but I shall 
give only the magnetisation 

^=Mo))U (10) 

for the d = 3,4 Ising model and the d = 3 Potts model (tables in this report. Further 

results and an analysis of the series will be published elsewhere or are available by email 
from the author. 
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1 
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5 
26 
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1071 
10010 
113636 
442220 
78367 
2424 
48 



2.5 50000 



Table 1: Number of connected diagrams for given dimension, maximum degree and number of 
spins. Diagrams with embedding number zero are not included. The last row gives the CPU time 
(p^ wall-clock time) needed. Note that the series with a smaller l^ax for a given dimension are 
not a prerequisite for the series with the largest Imax, but are listed for demonstration purposes. 
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Table 2: Magnetisation v — J^iViU^ for the d — 3 Ising model 
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Table 3: Magnetisation v — ViU^ for the d — A Ising model 
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Table 4: Magnetisation v = Y^i ^i^* for the d— 3 Potts model 
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